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摘 要 本 文 基于 离 焦 测 距 原 理 , 提出 一 种 利用 单 镜头 双 相 机 测量 运动 颗粒 三 维 速度 的 方法 . 入 射 光 经 分 束 立方 体 分 别 
进入 两 个 工业 相机 , 通过 调节 分 束 立方 体 与 相机 间 的 垫圈 厚度 获取 同一 对 象 两 幅 不 同 模糊 度 的 图 片 , 根据 其 模糊 差 实现 运 
动 颗粒 在 三 维 场 中 的 定位 , 再 结合 单 帧 多 曝光 方法 , 获取 其 三 维 速度 . 通过 拍摄 慢 速 搅拌 水 流 中 颗粒 运动 并 分 析 处 理 不 同 
运动 情况 的 颗粒 图 片 得 到 了 运动 颗粒 的 三 维 速度 ,经 计算 与 实验 考察 证 明了 其 可 行 性 . 该 方法 解决 了 离 焦 二 义 性 问题 与 速 
度 方向 二 义 性 问题 , 为 三 维 速度 的 测量 提供 了 新 的 方向 . 
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3-D Velocity Measurement of Moving Particles Based on Single-Lens 
Dual-Camera 
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(Shanghai Key Laboratory of Multiphase Flow and Heat Transfer in Power Engineering, School of Energy and Power Engineering, 
University of Shanghai for Science and Technology, Shanghai 200093, China) 
Abstract Based on the principle of depth from defocus, a kind of method for the measurement of 
three-dimensional (3-D) velocity of moving particle using the single-lens dual-camera was proposed. 
The incident light was dispersed into two industrial cameras through the beam splitter prism. Two 
images with different blur degree for the same object were acquired by adjusting the thickness of 
the washers between the prism and the cameras. The orientation of the moving particle in the 
three-dimensional field was obtained from the difference of blur degree between the images, and then 
the 3-D velocity of the particles can be calculated combined with the single-frame multi-exposure 
method. The images recording the movement of particles in the water flow with slow agitation were 
processed and analyzed which lead to the calculating results of the 3-D velocity field of moving parti- 
cles. The feasibility of the method proposed was verified through the comparison of the calculating 
results and the experimental results. The method solves the problems of the ambiguity of both de- 
focus and the direction of velocity, as well as indicating a new way for the measurement of 3-D velocity. 
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一 


0 引言 


多 相 流 在 燃烧 学 、 医 药学 等 各 个 学 科 领 域 中 有 
广泛 应 用 ,深入 了 人 解 多 相 流 机 理 有 赖 于 对 多 相 流 中 
颗粒 相 的 三 维 定量 测量 。 因 此 颗粒 场 的 三 维 测量 研 
究 具 有 重要 的 科学 意义 及 实用 价值 ， 如 药物 喷雾 的 
测量 、 循 环流 化 床 锅炉 气 固 两 相 流 流动 特性 的 研究 . 

目前 用 于 颗粒 场 三 维 速度 测量 的 技术 主要 有 : 
数字 全 息 技 术 (DH) 罩 、 激光 多 普 勤 测速 技术 


态 、 非 接触 测 基 等 优势 , 在 工业 领域 得 到 了 广泛 的 应 
用 . PIV 是 一 种 基于 杨 氏 条 纹 法 和 自 相 关 法 的 技术 ， 
其 可 在 极 短 的 曝光 时 间 内 连续 两 次 或 多 次 通过 CCD 
a ne 该 方法 在 二 维 速 
度 的 测量 上 已 相对 成 熟 ， 其 三 维 速度 主要 通过 与 全 

息 技 术 结合 , 即 数字 全 息 粒 子 测速 技术 (DHPIV) 加 ， 
或 者 采用 两 台 或 多 台 相 机 从 不 同 角 度 记 录 同 一 流动 


(LDV) 避 ee 束 技术 (PIV) 加 以 及 粒子 
跟踪 测 ; 速 技术 (PTV) 由 ,后 两 种 技术 由 于 其 全 场 、 瞬 
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区 域内 示 踪 粒 子 的 运动 情况 获取 加， 完成 颗粒 三 维 
速度 场 的 重建 . PTV 技术 从 原理 来 说 与 PIV 技术 近 
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似 , 相 比 于 PIV, 其 示 踪 粒子 浓度 较 低 , 避免 了 PIV 
技术 的 平均 效应 ; PTV 技术 可 对 单个 粒子 进行 跟踪 ， 
具有 更 高 的 准确 性 , 便于 由 二 维 测速 推广 到 三 维 测 
速 。 在 实际 应 用 中 每 种 方法 都 有 其 优势 和 局 限 ， 仍 
需 进一步 的 研究 探索 . 

20 世纪 90 年 代 初 迎 来 了 数字 图 像 处 理 技术 口 
的 大 发 展 。 该 技术 在 采集 图 像 过 程 中 通常 要 避免 产 
生 离 焦 ， 但 离 焦 时 产生 的 模糊 信息 能 反映 物体 的 深 
度 信息 ， 进 而 实现 三 维 位 置 和 速度 的 重建 。 因 此 本 
文 基 于 离 焦 测 距 原 理 ， 提 出 一 种 利用 单 镜头 双 相 机 
测量 运动 颗粒 三 维 速度 的 方法 .该 方法 解决 了 离 焦 
二 义 性 和 速度 二 义 性 问题 ， 操 作 人 简单 ， 且 可 行 性 得 
到 验证 ， 


1 实验 原理 


1.1 离 焦 测 距 原理 

物 方 远 心 镜头 加 离 焦 模糊 成 像 如 图 1 所 示 , 在 
像 方 焦 平 面 处 放置 一 个 直径 为 D 的 光 阑 , 使 物 距 发 
生变 化 时 物体 的 尺寸 大 小 可 以 保持 不 变 , 图 中 了 为 
物 点 , D 为 光 亲 直径 ，P' 为 像 点 . P 到 镜头 主 平面 
的 距离 为 物 距 uw， 光 闲 到 镜头 主 平面 的 距离 为 焦距 
f，P' 到 镜头 主 平面 的 距离 为 像 距 v。 


IS! IS, 岛 


二 
FN 


图 1 物 方 远 心 镜头 离 焦 模 糊 成 像 过 程 


Fig. 1 Defocus burring imaging process of object-space 


telecentric lens 


根据 高 斯 公式 二 + 了 = 于， 成 像 系统 的 焦距 / 
确定 后 , 物 距 wu 和 像 距 v 旦 一 一 对 应 关系 ， 只 有 相 
机 传感器 位 于 像 距 处 时 才能 得 到 最 清晰 的 像 ， 此 时 
为 最 佳 成 像 面 . 若 相机 传感器 偏离 最 佳 成 像 面 时 会 
产生 离 焦 模糊 ， 如 图 1 所 示 两 个 相机 传感器 分 别 位 
于 IS1 位 置 和 IS2 位 置 时 (s! > s?)， 此 时 传感器 偏 
离 了 最 佳 成 像 面 ， 会 产生 不 同 程度 的 离 焦 模 糊 ， 模 
糊 半 径 分 别 为 ri 和 72, ri 和 72 的 比值 根据 几何 关 

系 可 以 被 表示 为 
Ce 工 S2 干 UV 中 


结合 高 斯 公式 可 以 得 到 离 焦 图 像 的 物 距 计算 公 
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式 加， 
加 f(s1— Qs2) Qa>1,v<s2<s1 2) 
si—f—a(s—/) Q<1,s2<s1<v 

PE 


s1— f+a(ss—f) 


若 传 感 品 IS1 和 1S2 位 于 最 佳 成 像 面 的 同 侧 时 
选择 公式 (2); IS! 和 1S2 位 于 最 佳 成 像 面 的 异 侧 时 
选择 公式 (3). 公式 中 的 f, si 和 sz 可 以 通过 系统 标 
定 得 到 , a 通过 比较 图 像 的 离 焦 模糊 度 得 到 。 当 相 机 
传 感 铝 位 于 最 佳 成 像 面 时 会 产生 两 种 临界 情况 ， 即 
IS1 和 IS; 分 别 位 于 最 佳 成 像 面 时 ,由 此 可 以 推出 : 


= (8) 


其 中 , 51 为 IS2 恰好 位 于 最 佳 成 像 面 时 IS1 上 所 成 模 
糊 像 的 模糊 直径 ; 652 为 IS1 恰好 位 于 最 佳 成 像 面 时 
ISs 上 所 成 模糊 像 的 模糊 直径 . 若 2r1<01 且 2r2<0>， 
最 佳 成 像 面 位 于 双 相 机 传感器 之 间 ， 即 选择 式 (3)， 
车 271>01 或 2r?>02, 最 佳 成 像 面 位 于 双 相 机 传感器 
同 侧 ， 即 选择 式 (2). 
1.2 单 帧 多 曝光 方法 

为 获得 颗粒 场 三 维 速度 ,采用 单 帧 多 曝光 测量 
方法 , 如 图 2 所 示 . 第 一 个 外 部 触发 信号 输入 , 相机 
开始 接收 入 射 光 ， 进 行 第 一 次 上 曝光， 曝光 时 间 为 触 
发 信号 的 宽度 ， 经 过 一 段 时 间 后 ， 相 机 根据 设置 开 
始 进 行 下 一 次 上 曝光, 和 N 次 曝光 后 数据 集合 在 一 帧 图 
像 上 ， 形 成 单 帧 多 曝光 图 像 . 根据 实验 要 求 选择 合 
适 的 曝光 次 数 N、 曝 光 时 间 和 间隔 时 间 ， 大 曝光 次 
数 过 多 会 造成 噪声 增加 ， 上 曝光 次 数 过 少 则 无 法 准确 


时 时 间 上， 中 
相机 曝光 | : 
曝光 时 间 


数据 读 取 


多 曝光 数据 
图 2 单 帧 多 曝光 方法 


Fig. 2 Single-frame and multiple-exposure method 


图 3 为 一 幅 单 帧 多 曝光 图 像 ， 由 图 可 知 ， 一 个 
颗粒 经 过 四 次 曝光 后 在 一 帧 图 像 上 成 像 ， 直 观 反 映 
了 颗粒 的 运动 轨迹 ， 总 曝光 时 间 已 知 ， 只 要 选择 其 


201710.00076v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


1694 工程 热 物 理 学报 38 卷 


中 两 个 位 置 就 可 以 得 到 颗粒 的 三 维 速度 和 方向 。 为 
减少 误差 , 选择 位 置 1 和 4, 通过 程序 处 理 得 到 颗粒 
在 位 置 1 和 位 置 4 的 三 维 坐标 , 根据 下 面 公式 可 获 
得 颗粒 的 三 维 速度 : 


入 二 
T_T4— Xl (6) 
mt mt 
Ay yy4—Y 
一 一 7 
™y mt mt ( ) 
人 = 
v= (8) 
mt mt 


其 中 , m 为 镜头 的 放大 倍率 ; t 为 从 第 一 次 曝光 开始 
到 最 后 一 次 曝光 结束 的 总 时 间 ; w 和 ws 分 别 为 颗 
粒 处 于 位 置 1 和 位 置 4 时 的 物 距 ， 


图 3 单 帧 多 曝光 图 像 
Fig. 3 The image captured by single-frame and 


multiple-exposure imaging 


2 单 镑 头 双 相机 系统 


单 镜 头 双 相机 系统 由 两 个 CCD 工业 相机 ( 帧 率 
16 fbs、 分 状 率 1628x1236、 像 元 大 小 4.4 pmx4.4 
Wm)、 一 个 物 方 远 心 镜头 (1x) 和 一 个 分 束 立 方 体 组 
成 , 如 图 4 所 示 . 入 射 光 经 分 束 立 方 体 分 离 为 5:5 的 
透射 光 和 反射 光 , 分 别 进入 两 个 CCD 工业 相机 , 通 
过 调节 分 束 立 方 体 与 相机 之 间 的 垫圈 厚度 获取 了 同 
一 对 象 的 两 幅 不 同 模糊 度 图 片 ， 而 双 相 机 的 同步 运 
行 则 需要 使 用 信号 发 生 器 同步 触发 . 

镜头 是 由 多 片 透镜 组 成 的 透镜 组 ， 添 加 分 束 立 
方 体 后 可 近似 看 作 在 原 透 镜 组 上 添加 了 新 的 透镜 ， 
原 镜头 的 参数 被 改变 ， 故 镜头 与 分 束 立 方 体 的 组 合 
可 视 为 一 个 新 镜头 .该 组 合 的 焦距 与 其 他 参数 较 原 
镜头 可 能 发 生 改 变 ， 需 要 通过 实验 对 镜头 进行 重新 
标定 。 实验 选取 直径 为 1000 pm 的 圆 点 标定 板 作为 
拍摄 对 象 ,在 相机 和 分 光 棱 镜 之 间 从 0 开始 依次 添 
加 垫圈 (厚度 1 mm、2 mm 或 5 mm), 用 游标 卡尺 测 
得 分 束 立 方 体 到 相机 之 间 的 距离 ， 并 移动 圆 点 标定 


板 的 位 置 ， 找 到 该 像 距 下 对 应 的 最 清晰 像 的 工作 距 
离 (标定 板 到 镜头 前 端面 的 距离 )， 此 时 可 得 到 两 个 
参数 , 即 工作 距离 w 与 分 束 立 方 体 到 相机 之 间 的 距 
离 y, 由 高 斯 公式 得 到 公式 : 
0 
w+L1: vy+L2> ff 


其 中 ，Z1 表示 从 镜头 的 前 端面 到 其 物 主 平面 的 距 
离 ，L2 表示 从 镜头 的 后 端面 到 其 像 主 平面 的 距离 。 
根据 以 上 公式 拟 合 得 到 本 实验 镜头 和 分 束 立 方 体 组 
合 的 参数 如 表 1 所 示 . 


(9) 


三 维 颗粒 场 协 素 灯 


图 4 单 镜头 双 相 机 系统 示意 图 


Fig. 4 Schematic drawing of single-lens dual-camera system 


表 1 镜头 参数 
Table 1 The parameters of lens 
f/mm Li/mm L2/mm 
原 镜 头 47.76 29.14 78.01 
分 东 立 方 体 + 镜头 45.59 27.67 121.66 


在 系统 搭建 中 需 选 取 合 适 数 量 和 厚度 的 垫圈 ， 
使 同一 物体 在 双 相 机 上 的 成 像 模糊 度 不 同 ， 便 于 后 
续 程 序 处 理 得 到 其 模糊 度 之 比 。 系统 确定 后 ， 移 动 
圆 点 标定 板 , 分 别 记录 圆 点 在 双 相 机 上 成 清晰 像 时 
的 工作 距离 wt 和 wz， 再 分 别 加 上 Li 的 值 即 为 双 
相机 最 佳 成 像 面 所 对 应 的 物 距 ww， 代入 高 斯 公式 可 
以 得 到 像 距 w, 即 图 1 中 的 s1 和 s?。 经 计算 ，s1 为 
126.46 mm, s2 为 122.26 mm。 


3 实验 结果 


为 验证 系统 的 可 行 性 ， 对 水 槽 中 的 颗粒 场 进行 
拍摄 . 图 5 为 运动 颗粒 速度 测量 实验 装置 图 , 其 使 用 
粒 径 大 小 为 100 num、 密度 为 1.18 g/ems 的 示 踪 粒子 
PMMA( 聚 甲 基 丙烯 酸 甲 酯 ) 作为 流 场 内 运动 颗粒 的 
拍摄 对 象 . 实验 时 使 用 玻璃 棒 搅 动 水 槽 内 的 水 ,使 
其 低速 运动 ， 选 用 讽 素 灯 作 为 背景 光源 ， 采 用 背光 
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拍摄 方式 ， 通 过 信和 号 发 生 需 实现 双 相 机 的 同步 拍摄 
及 单 帧 多 曝光 功能 ， 曝 光 时 间 依 次 为 1 ms、4 ms、1 
ms 和 1 ms, 间隔 时 间 为 20 ms。 

图 6 为 运动 颗粒 在 双 相 机 中 成 像 的 几 种 典型 图 
片 及 处 理 结果 . 在 图 6(a) 中 , 沿 运动 方向 颗粒 在 两 
个 相机 中 成 像 的 模糊 度 变化 不 同 ,在 传感器 IS1 上 
模糊 度 逐 渐变 小 ， 而 在 传感器 IS2 上 模糊 度 逐 渐变 
大 ; 图 6(b) 中 运动 颗粒 的 模糊 度 变化 相同 , 均 逐 渐 
变 大 ; 图 6(c) 中 运动 颗粒 的 模糊 度 基本 不 变 . 由 于 
第 二 段 曝 光 时 间 较 长 ， 导 臻 颗粒 的 遮光 时 间 长 ， 除 
存在 一 定 长 度 运 动 轨迹 外 ， 相 比 其 它 颗粒 图 像 其 灰 
度 值 也 较 大 ， 通 过 该 长 轨迹 图 像 位 置 的 判断 可 以 获 
得 颗粒 的 运动 方向 。 为 避免 运动 模糊 带 来 的 复杂 计 
算 ， 应 选取 颗粒 位 于 短 曝光 时 段 的 图 像 进行 三 维 定 
位 处 理 计算 。 本 文通 过 对 颗粒 运动 轨迹 图 像 (椭圆 
形 框 中 ) 的 第 一 点 和 第 四 点 分 析 计 算 ， 得 到 zx、y、z 
三 个 分 量 的 速度 、 二 维 合 速度 和 三 维 合 速 度 ， 计 算 


结果 如 表 2 所 示 ， 其 中 z 分 量 以 向 左 运动 为 正 ，Y 
分 量 以 向 下 运动 为 正 ，z 分 量 以 向 镜头 方向 运动 
为 正 . 


Fig. 5 Experimental setup for velocity measurement of 
moving particles 


(0) 两 个 相机 中 沿 运 动 方 向 颗粒 图 像 模 糊 度 基本 不 变 
图 6 运动 颗粒 三 维 速度 测量 结果 


Fig. 6 Experimental results for velocity measurement of moving particles 
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表 2 运动 颗粒 的 三 维 速度 
Table 2 3-D velocity of moving particles (mmy/s) 


(mm/s) 


图 像 (a) 图 像 (b) 图 像 (c) 
vz 18.26 1.18 34.28 
vy 4.60 20.49 4.20 
vz 20.90 -0.73 0.00 
i 18.83 20.52 34.54 
es 28.13 20.54 34.54 


由 表 2 可 知 ,图 像 (a) 中 运动 颗粒 的 二 维 合 
度 vey 与 其 三 维 合 速度 we,wz 相差 较 大 ， 表 明 颗 粒 
相对 于 成 像 平 面 有 较 大 的 运动 ， 对 应 了 两 幅 图 片 运 
动 颗粒 图 像 模糊 度 变 化 不 同 的 情况 ; 图 像 (b) 和 (c) 
中 颗粒 的 vey 和 we,wz 基本 相同 ,表明 其 颗粒 在 z 
方向 上 基本 没有 运动 ， 对 应 了 两 幅 图 片 运动 颗粒 图 
像 模糊 度 变 化 一 致 的 情况 ; 图 像 (b) 和 (c) 分 别 对 应 
了 颗粒 垂直 运动 和 水 平 运动 。 实验 结果 证 明了 该 方 
法 和 系统 在 运动 颗粒 三 维 速度 测量 方面 的 可 行 性 。 


4 结 论 


本 文 基 于 离 焦 测 距 原理 和 单 帧 多 曝光 方法 , 利 
用 单 镜头 双 相 机 系统 测量 运动 颗粒 的 三 维 速度 ， 并 
简 述 了 高 焦 测 距 公 式 中 各 个 参数 的 获取 方法 ,通过 
拍摄 慢 速 搅拌 水 流 中 颗粒 的 运动 ， 对 不 同 运动 情况 
的 颗粒 图 片 进行 了 处 理 分 析 , 证 明了 该 方法 的 可 行 
性 与 有 效 性 , 操作 简单 、 成 本 低 , 实现 了 运动 颗粒 在 
无 参考 条 件 下 的 深度 定位 以 及 速度 重 构 ， 为 三 维 流 
场 测量 奠定 了 基础 。 
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